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The random-field Ising model shows extreme critical slowdown that has been described by acti¬ 
vated dynamic scaling: the characteristic time for the relaxation to equilibrium diverges exponen¬ 
tially with the correlation length, lnr ~ /T , with ip an a priori unknown barrier exponent. 
Through a nonperturbative functional renormalization group, we show that for spatial dimensions 
d less than a critical value Odr — 5.1, also associated with dimensional-reduction breakdown, ip = 0 
with 9 the temperature exponent near the zero-temperature fixed point that controls the critical 
behavior. For d > doR on the other hand, ip = 9 — 2A where 9 = 2 and A > 0 a new exponent. 

At the upper critical dimension d = 6, A = 1 so that ip = 0, and activated scaling gives way to 
conventional scaling. We give a physical interpretation of the results in terms of collective events in 
real space, avalanches and droplets. We also propose a way to check the two regimes by computer 
simulations of long-range 1-d systems. 

PACS numbers: 11.10.Hi, 75.40.Cx 

I. INTRODUCTION 

Activated dynamic scaling^ is a phenomenological de¬ 
scription of the extreme slowdown of dynamics observed 
in some disordered or glassy systems: systems in the 
presence of a quenched random field)^ spin glasses in 
their ordered phased pinned elastic manifoldsand 
possibly supercooled liquids as they approach their glass 
transition^ According to this scaling, the dynamics in¬ 
volve thermal activation over barriers that grow with the 
typical length scale l. leading to a characteristic time at 
the scale I behaving as lnr^ ~ ( E/T)(.^ with ip an a pri¬ 
ori unknown barrier exponent. Activated scaling leads 
to a broad distribution of relaxation times, which shows 
up in the time or frequency dependence of the response 
and correlation functions, and has also consequences for 
the nonequilibrium dynamics ji&l 

The random-field Ising model (RFIM)^ is one system 
whose dramatic critical slowing down is expected to be 
described by activated dynamic scaling^ Its critical 
point is controlled, in the renormalization group sense, by 
a zero-temperature fixed point at which the “dangerously 
irrelevant” renormalized temperature is characterized by 
an exponent 6 > 0. The dangerous irrelevancy leads to 
a breakdown of the hyperscaling relation between criti¬ 
cal exponents and to anomalous thermal fluctuations, all 
controlled by the exponent 0 and further rationalized at a 
physical level by the “droplet scenario” The simplest 
droplet assumption would be to set ip = 9. Actually, this 
equality has been found in the dynamics of a simpler dis¬ 
ordered system, an elastic manifold pinned in a random 
potentialj^iiS but in the case of the RFIM there has been 


no attempt to compute the barrier exponent ip. 

The functional renormalization group (FRG) is a tool 
of choice to provide a theoretical treatment beyond phe¬ 
nomenology and compute the barrier exponent ip. In 
its perturbative form, it has been successfully applied to 
the dynamics of the pinned elastic manifolds ££ For the 
RFIM, as was shown for the (static) equilibrium behav¬ 
ior, the FRG must be nonperturbative. 11-13 In this work 
we therefore generalize the nonperturbative FRG (NP- 
FRG) approach to describe the critical slowing down of 
the RFIM. 

We find that the critical slowing down of the RFIM is 
indeed of activated type with two different regimes as a 
function of spatial dimension d. For d less than a critical 
value doR — 5.1, also associated with the breakdown of 
the d — > d — 2 dimensional-reduction property ) 11-14 the 
barrier exponent coincides with the temperature expo¬ 
nent, ip = 9, as in elastic manifolds pinned in a random 
potential (see above). On the other hand, for d > doR, 
ip = 9 — 2A, where 9 = 2 (the dimensional-reduction 
value) and A > 0 is a new exponent that is computed 
within the NP-FRG. At the upper critical dimension 
d = 6, one finds A = 1 around the Gaussian fixed point, 
so that ip = 0, and activated scaling gives way to conven¬ 
tional scaling. We stress that in the range 6 > d > doR 
where the main critical exponents describing the static 
behavior coincide with the dimensional-reduction predic¬ 
tions, the critical dynamics is nonetheless activated and 
that this feature is completely missed by perturbation 
theory which instead predicts conventional dynamic scal¬ 
ing, t ~ £ z with 2 ~ 2 + 2u j 15 i 16 
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II. MODEL AND DYNAMICAL FIELD THEORY 

As we are interested in the long-time, collective behav¬ 
ior of the RFIM, a coarse-grained field theory provides 
an appropriate starting point. The relaxation dynamics 
of the scalar field p xt is thus described by a Langevin 
equation (for simplicity we consider the case of a non- 
conserved order parameter, known as model Ai£) 


After averaging over the thermal noises and the disor¬ 
der, one obtains 

Z[J a j Ja\ = 

Dip a X>ip a e~ Sdvn ^ a ’‘ Pa ^ + ^ a + 

v a 

(4) 

where the (bare) dynamical action is 


dtPxt = -&b 




+ T) xti 


(i) 


where r] x t is a Gaussian random noise term with zero 
mean and variance ( r] xt ri x 't>) = ZTVLbS^ (x — x')5(t — t'). 
The ’’action” or effective Hamiltonian S'fiyj] is given by 


h + J]= S B [<p]- / [h{x) + J x \p x , 

J X 

= J | \ldxWx] 2 + 


( 2 ) 


where f r = f d d x, J x is an external source, h x is a ran¬ 
dom “source” (a random magnetic field) taken with a 
Gaussian distribution characterized by a zero mean and 
a variance h x h x > = AbS^(x — x'). 

The generating functional of the multi-point and multi¬ 
time correlation and response functions can be built as 
usual by following the MSR formalism] 18 i 19 After intro¬ 
ducing an auxiliary “response” field ( p xt and taking into 
account the fact that the solution of Eq. m is unique *22 
one obtains the “partition function” 


Zh,r)[J, J] 


J VtpVip exp { 



dtPxt + 


Spxt 


f]xt hx 


(JxtPxt T 


( 3 ) 

where we have used the Ito prescription (which amounts 
to setting to 1 the Jacobian of the transformation be¬ 
tween the thermal noise and the field), 20 

The conventional route for studying the dynamics of 
disordered systems is then to average the partition func¬ 
tion in Eq. m over both the thermal noise and the 
disorder and to take advantage of the property that 
Zh,r/[j = 0, J] = 1.— However, in previous NP-FRG 


work on the RFIM*** — it was shown that the key point 
for taking relevant events such as avalanches and droplets 
into account is to describe the full functional dependence 
of the cumulants of the renormalized disorder, a point 
that is overlooked by the standard replica, superfield or 
dynamic formalisms. The most convenient procedure to 
obtain this full functional dependence is to introduce 
copies or replicas of the system: the copies have the 
same disorder h but are coupled to distinct sources, in 
contrast with the usual replica trick . 11 ’ 13 We therefore 
combine dynamics and replicas/copies. The latter are 
now characterized not only by distinct sources but also 
by independent thermal noises* 22 


Rdyn\{Pa->Pa\\ — ^ ] / Pa,xt\S)tPa,xt T(p a ^ a 
a Jxt ^ 

*[<Pa] \ / * 

} o L ^ a > 

a,xt J z , J xtt' 


SS B [p, 

Sip, 


( 5 ) 


xt^Pb,xt' 


and where we have set fls = 1; In Z is the sought generat¬ 
ing functional of the response and correlation functions. 

In the long-time limit, the relaxation toward equilib¬ 
rium satisfies, in addition to the causality requirement, 
an invariance under time translation (TTI) and a time- 
reversal symmetry (TRS). 20 The latter in turn implies 
the fluctuation-dissipation theorem3^i 2 4 The TRS cor¬ 
responds to an invariance of the theory under the si¬ 
multaneous transformations t —> —t, p a -A p a , and 
Pa t ip a (1 /T)d t pa~ 


III. NONPERTURBATIVE FUNCTIONAL 
RENORMALIZATION GROUP 

The theoretical formalism we use to describe the long¬ 
time, long-distance physics of the RFIM near its critical 
point is the NP-FRG. We have generalized the formal¬ 
ism developed for the (static) equilibrium properties of 
the RFIMii - — 3 by combining it with the approach put 
forward by Canet et al. 2 - for the critical dynamics of the 
Ising model in the absence of quenched disorder. 

To apply the NP-FRG formalism to the above dynam¬ 
ical field theory, we introduce an infrared (IR) regulator 
A Sk to the action (0, whose role is to suppress the inte¬ 
gration over slow modes associated with momenta |g| < k 
in the functional integral: 11 ’ 13 ' 25 ’ 26 

AS’ fe [{$a}] = ^ / tr\^2<f> a}Xt R k (\x-x'\,t-t')$T x , t , 

^ J xx'tt' L a 

+ ^ ^ > ^a,xt^k{\ x — x \, t — t )^b,x't' » 
ab 

( 6 ) 

where $ a = ( p a ,p a ), its transpose, and Rfc and 

Rfc are symmetric 2x2 matrices of mass-like IR cutoff 
functions that enforce the decoupling between fast (high- 
momentum) and slow (low-momentum) modes in the par¬ 
tition function. Following Ref. [27] it proves sufficient 
to control the contribution of the fluctuations through 
their momentum dependence and take Rfc,n = -Rfc ,22 = 0, 
Rfcq 2 = Rfc^i = Rk{x—x'), and Rk,n =^Rk,12 = Rfc, 21 = 

0, Rfc, 22 = Rfc( x—x') where Rk(q 2 ) and Rfc(g 2 ) are chosen 
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(in Fourier space) such that the integration over inodes 
with momentum |g| < k is suppressedTo avoid 
an explicit breaking of the underlying super-rotations of 
the theory^ we takei^ 

n /„2'\ „ dR k (q 2 ) /r7 , 

Pk(q ) oc QqZ ' v) 

Note that the above choice of IR regulator satisfies the 
TRS, a crucial property. 

Through this addition Z[{J a }] is replaced by a k- 
dependent quantity, Z k [{J a }\, where J a denotes ( J ai J a ). 
The central quantity of the NP-FRG is the “effective av¬ 
erage action” which is the generating functional 

of the 1 -particle irreducible correlation functions at the 
scale k. It is defined (modulo the subtraction of a regula¬ 
tor contribution) from In Z k \_{J r 0\\ via a Legendre trans¬ 
form: 

r fc [{$ a }]+in z k [{j a }] = J2 [ teJa,xt$Z xt -&s k [{t> a }], 

a Jxt 

( 8 ) 

where $ a = (0 a ,0 a ) now denotes the “classical” (or av¬ 
erage) fields with 0 a , xt = 6 In Z k /6 J a , xt = (y a ,xt) and 
0a, Xt = 5 In Z k /5J a ,xt = (0a,xt) j the trace operation is 
over the 2 components of <h a and J a - 

Expansions in generalized cumulants are then gener¬ 
ated by expanding the functionals in increasing number 
of unrestricted sums over copies, 

00 i—iiP- 1 

r fe [{<U] = E L Er- E r fep [d> ai ,•■•,<&«,] (9) 

p— 1 ai---a p 

where V kp can be formally expressed as 

r kp = [ 0a 1 ,x ltl ---0a P ,x p t p lk P -,x 1 t 1 ,-,x p t p (10) 

with 7 kp a functional of the fields 'I’ai.tn • • • , &a p ,t p and 
of their time derivatives, $ ai , tl , ■ ■ ■ ,dt p ®a p ,t p , 9 > 1 - 
When the fields are chosen uniform in time with <j> a ,xt = 
0 a , x and 0 a ,xt = 0 , the 7 fc p ’s reduce to the cumulants 
of the renormalized random field at equilibrium already 
studied in Refs. pTl4l3l |. For generic fields, the additional 
contributions then represent kinetic terms£ 

The functional T k satisfies an exact RG equation 
(ERGE) that describes its evolution with the IR cutoff 
k^. 

W{$ a }] = iTr{(^R k )(rf [{dTJl + Rk)- 1 }, (11) 

where the trace is over space-time coordinates, copy in- 
dices and components, and F).' is the the matrix formed 
by the second functional derivatives of T k . (In what fol¬ 
lows, superscripts within a parenthesis are used to in¬ 
dicate derivatives with respect to the appropriate argu¬ 
ments.) By inserting the expansion in increasing num¬ 
ber of sums over copies and proceeding to the associated 
algebraic manipulations, one then derives an infinite hi¬ 
erarchy of ERGE’s for the generalized cumulants T kp or, 
alternatively, for the functionals ^ kp . 


IV. NONPERTURBATIVE APPROXIMATION 
SCHEME 


One cannot hope to solve the infinite hierarchy of func¬ 
tional flow equations for the 'y kp , s exactly, but one can de¬ 
scribe the long-distance physics of the problem by means 
of a nonperturbative approximation scheme. We com¬ 
bine the minimal truncation of the effective average ac¬ 
tion already shown to successfully describe the equilib¬ 
rium critical behavior of the RFIMi^ with an account of 
the dynamics through a truncation of the expansion in 
kinetic coefficients that allows us to describe the charac¬ 
teristic relaxation time. By taking into account the TRS, 
we arrive at the following ansatz: 


7fci;.ri [*J1] — , U k (0 x t ) V ~Z k (0 x 0)(d x 0 x t)~ 

0xt L Z 

“b X k (0 x ±)(dt0xt T0xt) 


( 12 ) 


7fe2;xiti,S2t2 [^1; ^ 2 ] — S^ d \xi — X 2 ) A k (0i, Xl t 1 , 02,x x t 2 ) 

(13) 

while the 7 fe p ’s with p > 3 are set to zero. Note that 
this ansatz describes the characteristic relaxation time 
but not its distribution: to do this the next orders of the 
truncation of the expansion in kinetic coefficients would 
be required, as discussed in Refs. [H0|- 

The next step is to derive the RG flow equations for 
the functions contained in the ansatz from the ERGE’s 
for the 7 /cp’s. As already mentioned, we work in the Ito 
discretization scheme and the corresponding prescription 
can be systematically implemented in the NP-FRG equa¬ 
tions by following the simple procedure developed in Ref. 
[27}. The derivation is tedious but straightforward and 
more details are given below and in the appendix. The 
output is a set of coupled flow equations for 3 functions 
of one field, U' k (0), Z k (0), X k (0), and one function of 2 
fields A k (0i,02). As a result of the TRS, the renormal¬ 
ized kinetic function X k does not enter the flow of the 
static ones (U' k , Z k and A k ). 

The flow equations involve the renormalized propaga¬ 
tors at scale k evaluated at the lowest order of the expan¬ 
sions is sums over copies and for fields that are uniform in 
space and time, $> a ,xt = (0a, 0). In Fourier (momentum) 
space, these propagators are expressed as P k ,ab(q', t, t') = 
P k (q-, 0 a ;t,t')5 ab + P k (q ; 0 a ,0b\t,t') where the 2x2 ma¬ 
trix Pfc has a structure following from causality, TTI and 
fluctuation-dissipation theorem, with 

P k 2 (q; 0 ; t' -t) = 0(f' - t)X k (0)~ 1 e~ (14) 

the response function, P 21 (t' — t) = P k 2 (t ~ f) , P^ 1 the 
2-time disorder-connected correlation function given by 

?l\q ] 0 ] t'-t) = T[Z k (0)q 2 +R k (q 2 ) + U'a0)]- 1 e -&T 

(15) 
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and P% 2 = 0; the characteristic relaxation time is defined 
as 


Tk(q-, <t>) 


_ x k (<t>) _ 

[Z fc (0)g 2 +iW)+ !/"(</>)] ' 


(16) 


In addition, the only nonzero component of Pfc is 


H 1 (q> <i>u<h'A -t) = 

fPl 2 (g; ^; t) f P 21 (g; </> 2 ; t') [A* (</>i, </> 2 ) - R k (q 2 )] 

(17) 

which, in this truncation, is simply the static (equilib¬ 
rium) disorder-disconnected correlation function. 

Finally, to study the vicinity of the relevant zero- 
temperature critical fixed point 2 the NP-FRG equations 
are cast in a dimensionless form by introducing appro¬ 
priate scaling dimensions: (f) ~ fc( d - 4 +W 2 j z k ~ k~ v , 
U' k ~ k(. d ~ 2 v+v)/ 2 , A/. ~ k~( 2ri ~^ and the renormalized 
temperature T k ~ Tk 6 , where 9 = 2 + r] — rj > 0 . 11 ' 13 We 
express the results in terms of the dimensionless fields 
ip = Vl + V2 and Sip = V2 ~ Vl . With lowercase letters de¬ 
noting dimensionless quantities, one can formally write 
the flow equations for the static quantities as 

kd k u' k (ip) = /3„' 0 {ip) + T k j3 u 'i(<p ), 

kd k z k {ip ) = Pzo(q>) + T k fi z i(tp ), (18) 

kd k 6 k (ip,Sip ) = f3so{<P,Sip) +T k /3si(<p,Sip ), 


where the beta functions depend on u’ k , z k , S k , their 
derivatives, and on the (dimensionless) cutoff functions. 

These equations generalize to nonzero temperature 
those given in Ref. [Hi: /3 u > 0 (<p), /3 2 o(t), and pgo(ip,Sip) 
coincide with the zero-temperature beta functions explic¬ 
itly given in this reference. The beta functions /3 u ’i{ip) 
and P z i(ip) are regular functions whose expression is unil¬ 
luminating and is not given here. Finally, /3gi(ip,Sip) is 
given by 


Psi(<P,8<p) = ~fd a r{q 2 ){\p k (q;<p + S<p) 2 + sym] x 

[4 02 VM + ^ 0 \ip,Sip)\ + 2[p k (q; ip + Sip) 2 - sym] x 

S[ U \v,Sip)+ 2S^ 1) (ip,Sip)-ff-\p k (q;ip +Sip) 2 + sym] 

+ 26 { k °\ip, 5ip)-^-\p k (q-, ip + Sip) 2 - sym] 

^ (19) 

where q = q/k, f~ = f d d q/(2n) d , the dimensionless cut¬ 
off function f{q 2 ) is defined through R k (q 2 ) = Z k q 2 r(q 2 ) 
and d s r(q 2 ) = —[r]q 2 f{q 2 ) + 2q 4 r'(q 2 )\ is a symbolic no¬ 
tation for the term obtained from kd k R k (q 2 ). Simi¬ 
larly, one defines f(q 2 ) from R k {q 2 ) = A k r(q 2 ) but it 
is simply related to f via r(q 2 ) = —8^. [q 2 r(q 2 )] from 
Eq. 0. The dimensionless hat propagator is given by 
Pk{q;<p) = ( q 2 [z k (ip) + r(q 2 )\ + u'l{ip))~ 1 . Finally, sym 
denotes a term obtained by changing Sip in —Sip. 


When T = 0, it was previously found that the fixed- 
point solution displays two regimes— 

(1) for d < dDR — 5.1, a “cusp” in \Sip\ is present in 
the fixed-point function 5* when Sip —» 0: 

S*(ip, Sip) = S*(ip, 0) - 5* t a(ip)\5ip\ + 0(5ip 2 ) (20) 

with ^ 0. This cusp is associated with the presence 
of avalanches on all scales at the critical point— 

(2) For d > dDR the fixed-point function <5* is “cus¬ 
pless” , which ensures that the d —> d — 2 dimensional- 
reduction property of the (static) critical exponents^ is 
valid (and that the super-rotation is not spontaneously 
broken along the RG flow). Avalanches are still present 
but lead to only a subdominant cusp^ 

S k (q>, Sip) = S*{ip, 0) - 5k,a{v)\Sip\ + 0(5ip 2 ) (21) 

when k —> 0, where S kja goes to zero as 

6kA<f>) ~ ( 22 ) 

with A > 0 characterizing the (diverging) number of 
spanning avalanches— AL 


V. THERMAL BOUNDARY LAYER 

Describing the critical slowing down requires a nonzero 
temperature and additional care is needed. The beta 
function for S k shows a nonuniform convergence when 
k —> 0 and Sip —> 0 and for a nonzero T the cusp is 
rounded near Sip = 0 in a “boundary layer” of width 
Sip ~ Tfc (see also Refs. [36!] for the case of the elastic 
manifold in a random environment). 

For T ^ 0 the beta function of m the limit y 0 
and Tk —> 0 can indeed be written as 

7.0 x t..\ , a lk(<p) 8 2 [Sk{ip,Sip) - 4(t,0)] 2 

kd k S k {ip, Sip) —Ps,reg{ip) P ^ d(Slf) 2 

- Tfca 2 fe(v 5 )4° 2) (t> Sip) 

(23) 

where Ps,reg is the contribution that is independent of 
the derivatives of S k with respect to Sip; a\ k and a 2 fc 
are regular functions obtained from the static functions: 
ai kiip) is the prefactor of the anomalous contribution in 
f3go(ip,5ip) whose limit when —>• 0 is nonzero only in 

the presence of a cusp, 

aik{<p) = \f d sf{q 2 )p k {q\ ipf , (24) 

and a, 2 k{ip) is the prefactor of the potentially singular 
piece in Psi(ip,Sip), 


a 2k (<p) = \f d sr(q 2 )p k {q; ip) 2 . (25) 
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The cusp in \8<p\ that is present in T = 0 is rounded at 
finite temperature because of the last term in Eq. (1?51) . 
Instead, 8k develops a “thermal boundary layer”, 

6k(<P, M = 4(<A 0) + T k b k (ip , y=7f-) + 0(T 2 , Sip 2 ) 

(26) 

when Xfc, 8ip —> 0. It is easy to derive that the solution 
has the explicit form 


The solution of Eq. (l30l) when k —> 0 is then of the 
form 

F k(‘P) = 7f-ek + -L=9k(‘p) + 0{ 1) (31) 

with efc independent of ip. We find that e k —> e* > 0 for 
d < dun and e k — > 0 as k 2X for d > duR [compare with 
the term in 1/Tfc in Eq. (l30l) ]. 


fo(<P,V ) 


a2*(p) 

au(<p) 



ai*{ip) 2 8k,a{p) 2 

a 2 *(ip) 2 



(27) 


where a p * are the (nonzero) fixed-point functions and 
8k,a{p) behaves differently when k —> 0 for d < dpR and 
d > dDR (see above). 


VI. ACTIVATED DYNAMIC SCALING 


We now turn to the NP-FRG equation for the kinetic 
term X k ((f>). It is obtained from the renormalization pre¬ 
scription 


X k (<f>) = 


rp „ 2 TfcljSCt (</*) </>) 

J- dtp 


0=0 


(28) 


and is given in graphical terms in the appendix. After 
having introduced the dimensionless quantities, it can be 
rewritten as 


kd k X k (ip) = /?xoM + T k /3 X i{ip) (29) 

where /3xo(p) is given in the appendix and f3xi(p) is a 
regular function that leads only to subdominant terms 
near the fixed point. 

Note that we have kept X k itself in a dimensionful 
form. In the case of a conventional critical slowing 
down, 17 one introduces a dynamical exponent z such that 
the characteristic relaxation time scales as r k ~ k~ z near 
the fixed point. The kinetic term then has dimension 
X k ~ k~ z+ ( 2 ~ r, \— However, in the presence of a nonzero 
random-field strength, where one anticipates an uncon¬ 
ventional activated dynamic scaling, one should rather 
focus on F k = ln.Yfc (where if needed X k can be made 
dimensionless inside the logarithm by dividing by a k- 
independent factor). 

By inserting the boundary layer solution in Eqs. m 
and (IA3I) (see the appendix) and working at the dominant 
orders when T k — > 0 (and k —> 0), it is easy to derive the 
flow of F k {ip) = \nXk{4>), which reads 

kd k F k {ip) = - — ^^-ipF^ip) + c lk (ip) + c 2k (ip)F^<p)+ 

C 3 k{p)[F k (ip) + F k {ip) 2 } - -L + o{T k ), 

J-k CL2k[P) 

( 30) 

where the c pk s are regular functions of ip whose expres¬ 
sions are given in the appendix. 


VII. RESULTS 

To support the analytical solutions provided above and 
compute the exponents we have numerically solved the 
NP-FRG equations for a wide range of dimensions be¬ 
tween 3 and 6. To do so, we have discretized the fields 
on a grid and used a variation of the Newton-Raphson 
method. For the cutoff function f k {q 2 ) we have chosen 
the same form as in previous work and optimized the 
parameters by stability considerations^ 3 

First, the numerical solution confirms the behavior of 
F k (ip) given in Eq. (1311) . We illustrate this point by 
showing the flow of v / 7fc[-Ffc(y>) — Ffc(O)], which should 
asymptotically converge to gk(p) — g k { 0), for the case 
d = 4.4. We can see from Fig. [I] that the fixed-point 
function is indeed well-behaved. 



FIG. 1: Evolution of y/Tk[F k (ip) — F k (ip o)] for RG times s = 
lnfc from 0 (initial condition where the function is chosen 
equal to zero: red line) to -7.1 (essentially, the fixed point: 
blue line). Here ipo = 0 (but other choices lead to the same 
asymptotic function). The function is illustrated for the case 
d = 4.4 and a bare temperature T = 0.1; the minimum of the 
effective potential is then for p ~ ±0.085. 


In Fig. [2] we further illustrate the dominant 1/T*, de¬ 
pendence of F k for the two cases discussed above with 
F k ~ 1/Tfc for d = 4.4 < duR and F k ~ k 2X /T k for 
d = 5.4 > duR- 

As a result of Eqs. and ED: the asymptotic be¬ 
havior of lnTfc goes as e k /T k with e k given above and 
Tfc ~ Tk®, so that the characteristic relaxation time r k 
follows the activated dynamic scaling with 


Tfc 


e T 


(32) 
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FIG. 2: The flow of \nFk(ipa) with the RG “time” s = ln(fc/A) 
with the UV scale A = 1; tpo is (arbitrarily) chosen near the 
minimum of Uk{ip), for d = 5.4 and d = 4.4. The thin full 
lines denote the expected asymptotic slopes, ip = 8 ~ 1.909 
for d = 4.4 and ip = 8 - 2A « 0.855 for d = 5.4. 



FIG. 3: The barrier exponent ip (full black line) and the tem¬ 
perature exponent 8 (dashed red line) versus dimension d; 
8 ~ 1.49 in d = 3 and 1.84 in d = 4, in excellent agreement 
with simulation results— i 


for its leading behavior when k —> 0 (we have dropped the 
prefactor and subdominant terms in the exponential). 

From the above, we find ip = 6 for d < doR and ip = 
6 — 2A with 9 = 2 (due to dimensional reduction) and 
A > 0 for d > don- At the upper critical dimension d = 
6, one finds that A = 1 around the Gaussian fixed point 
so that ip = 0— Activated dynamic scaling thus gives 
way to conventional dynamic scaling/^ a scaling which 
is for instance found in the mean-field, fully connected, 
version of the model. The barrier exponent 'ip is shown 
as a function of d in Fig. [3] 


VIII. DISCUSSION 

The above results can be physically interpreted by in¬ 
voking the collective events that are present at critical¬ 
ity. At T = 0 these are avalanches, which have a frac¬ 
tal dimension df = (d + 4 — fj)/2 — 2A with A = 0 for 
d < doR and A > 0 for d > doR^ For T > 0 but very 


small, there emerge from these avalanches, which cor¬ 
respond to an exceptional degeneracy between ground 
states, low-energy excitations corresponding to quasi- 
degeneracy, known as droplets. The result for the bar¬ 
rier exponent can then be rationalized by assuming that 
with a probability T/L ( 0-2A ) there are critical samples of 
size L with such a quasi-degeneracy and that the whole 
energy landscape of the quasi-degenerate system has a 
unique scale, then given by L^ e ~ 2X \ 

The droplet picture also predicts anomalous static 
thermal fluctuations of the field (or magnetization) — 
In the rare samples, the magnitude of the magnetiza¬ 
tion fluctuations goes as L d f , so that the pth cumulants 
associated with the thermal fluctuations of the magneti¬ 
zation, [L~ d < (f ip x — < f <p x >) 2 >] p , h as an anoma¬ 
lous scaling, oc (T / L 8 ~ 2 X )L^ 4 ~ P ~ 2X ' >P . We have checked 
the validity of the scaling for p = 2 from the NP-FRG 
equations, along the lines detailed in Ref. [l^ |. 

Finally, we conclude by proposing a way to directly 
check the two regimes, ip = 9 and ip = 9 — 2A with A > 0, 
in computer simulations. As recently shown^ the 1 -d 
RFIM with long-range power-law interactions oc |*| d+<T 
has a critical value a c ~ 0.379 around which the change 
of regime should be observed. In the mean-field region, 
for a < 1/3, one should also find that the relaxation is 
no longer activated but follows conventional scaling with 
the mean-field dynamic exponent z = o. 


Appendix A: NP-FRG flow for the kinetic term 

The flow of the kinetic term Xk(<p ) is obtained from the 
renormalization prescription in Eq. (1281) and from the 
ERGE for 7 fci ;x t. In graphical terms the flow equation 
reads 



(Al) 


where a cross in the circle denotes a vertex Xk{4>), the 
lines denote static propagators Pk and the dotted lines 
the disorder vertex Ak((pi, $ 2 ) (after having taken the 
needed derivatives, one sets </>i = (p 2 = &fc is a short¬ 
hand notation to indicate a derivative acting only on the 
cutoff functions, i.e., dk = dkRkd/SRk + dkRkd/5Rk- For 
implementing Ito prescription we have followed the trick 
devised by Canet et al.—, which amounts to shifting the 
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time dependence of the response field by an infinitesimal 
amount in the renormalized response functions. 

With the help of the dimensionless quantities, Eq. 
ED can be rewritten as 

kd k X k {ip) = @xo(<p) + T k p X i(<p) (A2) 


where 


C 2 k{p) = 2 


/ Pk{q;<pY 

Jq 


d s f{q 2 ){q 2 z k {ip)+u^ ) {p)) 


+ d s f(q 2 )^-6 ( j} 0 \ip,0)+3p k (q;ip)[5 k (ip,0)~r(q 2 )] 


x [Q 2 z' k (ip)+u[ ■ 


0(T k ). 


(A5) 


Pxo(<fi) = - — Y^<pX' k {tp) + i 


- 2d s f{q 2 )p k {q;p)‘ 2 


x A k {p) - 4 p k (q; pf [q 2 z k (<p) + wfV)] ( -2 [d s r{q 2 ) 
+ 3 d s r{q 2 )p k {q; ip)(5 k (ip, 0) - r(q 2 ))] X k (<p) - 3 d s f{q 2 ) 
x p fc (g; ^)4 1,0) (v, 0)X k (<p)) - 4 p k (q-, p) 4 X k {p) [d s f (q 2 ) 


C3k = — 


+ 4d s r{q 2 ){5 k {p, 0) - f{q 2 ))p k {q\ ip)\ [q 2 z k {p) + u[ ; (p)] 

- 2 d s r{q 2 )p k {q; pf [2((5(<p, 0) - r{q 2 ))X k {p) + 4 6^'°\p, 0) 


d s f{q 2 ) + 2d s r{q 2 )[5 k {p, 0) - r{q 2 )]p k {q; ip 

Pk(q ; v) 2 ~T k d s r{q 2 )p k {q ; ip) 2 . 

J q 


(A6) 


x' k {p) + (- 4 0,2) (<p> o) + 4 2,0 V o ))x k (<p) 

(A3) 

where d s f{q 2 ) is a symbolic notation for [2rj — rj\[r(q 2 ) + 
g 2 r'(g 2 )] + 2[2g 2 r'(<j 2 ) + g 4 r"(g 2 )]. The regular term 
0x\ijp) leads only to subdominant terms near the fixed 
point and is not given here. 

The flow of F k {ip) = In X k (</>) at dominant orders when 
Tfc —> 0 is given in Eq. (15(1 and the (regular) functions 
c pk {p) are expressed at the relevant order in T k as 


ci k(<f) = 7 } J [ ~ 2 d s r(q 2 )p k {q; p) 4 [q 2 z' k {p) + {p)] 

+ d s f{q 2 )p k {q\ ip) 3 ( - S (20) {p, 0) + 2 p k {q- p)[q 2 z k {p) 
+ u< k \ ( P )] [ 3 4 10) (^.°) - (4(^,o) - f(q 2 )) 


x W 2 4M + ' 


+ 0{T k ), 


(A4) 
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